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The  ONR  Sea  Loads-Ship  Response  (SLSR)  program  includes  research  areas  related  to  nonlinear 
hydrodynamics,  nonlinear  dynamics,  structural  fatigue,  elastic  and  plastic  structural  deformation, 
and  a  probabilistic  or  reliability-based  analysis  of  ship  structural  design.  Due  to  the  diverse  and 
multi-disciplinary  nature  of  the  project,  program  researchers  were  brought  together  at  the 
University  of  Michigan  to  discuss  the  direction  of  their  current  and  future  research;  the  goal 
being  to  achieve  a  high  level  of  coordination  between  the  various  efforts.  This  report  contains 
copies  of  the  presentations  made  at  the  workshop. 

Thirty-one  participants  from  various  academic  institutions,  government  laboratories  and  offices, 
and  commercial  companies  attended.  Presentations  representing  the  state-of-the-art  were  made 
in  the  areas  of  hydrodynamic  loading,  structural  analysis,  design  reliability,  and  simulation-based 
design. 

Experts  in  hydrodynamics  (SAIC,  MIT,  AMI,  and  UofM)  explained  that  by  using  various 
nonlinear  or  partially  nonlinear  models,  computer  codes  are  capable  of  determining 
hydrodynamic  loads,  excluding  bottom  impact  or  flare  slamming,  in  random  seas.  However, 
given  the  success  of  recent  planing  hull  studies,  the  extension  of  planing  hull  hydrodynamics  to 
the  impact  problem  should  be  straightforward  thus  allowing  for  the  complete  hydrodynamic 
loads  time  history  in  extreme  seas  to  be  made.  From  these  time  histories,  the  design 
hydrodynamic  and  inertial  loading  events  can  be  determined. 

Structural  experts  (CDNSWC,  NAVESEA,  Ross  and  McNatt,  and  ABS)  explained  how  the 
hydrodynamic  and  inertial  loads  are  currently  estimated  and  used  in  the  structural  design  of 
ships,  both  naval  and  commercial.  Due  to  the  complexity  of  a  ship's  structure  and  the  need  for 
timely  engineering  answers,  hydrodynamic  load  modeling  is  generally  simpler  than  that  available 
as  described  by  the  previous  hydrodynamic  experts.  It  was  agreed  that  hydrodynamic  and 
structural  analysis  code  integration  is  a  high  priority  of  the  SLSR  project  and  means  for 
achieving  this  integration  were  identified. 

Finally,  experts  in  reliability,  virtual  reality,  and  simulation  (UofC,  NRC,  and  UofM)  gave 
examples  of  how  the  product  of  the  SLSR  program  could  fit  into  a  larger  computer  environment 
where  simulation-based  designs  incorporating  probabilistic  methods  would  be  possible. 

In  summary,  the  workshop  was  one  of  the  few  times  where  researchers  of  the  disparate 
disciplines  were  brought  together  to  develop  a  coordinated  program  for  ship  structural  design. 
Through  the  spirited  discussions,  a  new  awareness  of  the  problems  facing  the  different  fields  was 
formed  and  in  this  respect,  the  workshop  must  be  considered  a  success. 
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Workshop  Introduction 
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8:30  -  8:45  Welcome 

Prof.  Michael  G.  Parsons,  Naval  Architecture  and  Marine  Engineering, 
Associate  Dean,  College  of  Engineering,  University  of  Michigan 


8:45  -  9:00  Workshop  Focus 

Dr.  Peter  Majumdar,  Office  of  Naval  Research 


9:00  -9:15  NAVSEA  Initiative  in  Wave  Loads  Predictions 

Mr.  Allen  H.  Engle,  Naval  Sea  Systems  Command 

Hydrodynamics 

9:15  -  9:50  Large-Amplitude  Motion  and  Wave-Load  Predictions  for  Ship  Design 
Assessment 

Dr.  Nils  Salvesen,  SAIC 


9:50  -  10:25  Nonlinear  Ship  Motions 

Prof.  Paul  D.  Sclavounos,  Ocean  Engineering,  Massachusetts  Institute  of 
Technology 

10:25  -  10:40  Break 


10:40  -  11:15  Prediction  of  Nonlinear  Loading  of  Flared  Bodies  Using  a  Numerical 
Towing  Tank 

Dr.  Brian  Maskew,  Analytical  Methods,  Inc. 


11:15-11:50  Fully  Nonlinear  Hydrodynamic  Loads  Using  De-Singularized  Methods 

Prof.  Robert  F.  Beck,  Naval  Architecture  and  Marine  Engineering, 
University  of  Michigan 

1 1 :50  -  12:25  Loads  Associated  With  the  Hydrodynamic  Impact  of  Flat  Wedges 

Prof.  William  S.  Vorus,  Naval  Architecture  and  Marine  Engineering, 
University  of  Michigan 


12:25-1:25 


Lunch 


1 :25  -  2:00  Non  Hydrodynamic  Forces  on  High  Speed  Vessels 

Pro*  .  vanin  W.  Troesch,  Naval  Architecture  and  Marine  Engineering, 
University  of  Michigan 

Structures  and  Design 

2:00  -  2:35  Ship  Structures  and  NAVSEA 

Mr.  Jerome  P.  Sikora,  CDNSWC 

2:35-3:10  Integrated  Ship  Structural  Design  Methodology 

Mr.  Tobin  R.  McNatt,  Ross  and  McNatt  and  Prof.  Owen  Hughes, 
Aerospace  and  Ocean  Engineering,  Virginia  Polytechnic  Institute  &  State 
University 

3:10-3:25  Break 

3:25  -  4:00  Probabilistic  Loading  of  Shi^  Stn.ctures  by  Slamming 

Prof.  William  Webster,  Naval  Architecture  and  Offshore  Engineering, 
University  of  California,  Berkeley  (for  Prof.  Alaa  Mansour) 

4:00  -  4:35  Dynamic  Loading  Approach  for  Analyzing  the  Ship  Structure 

Dr.  Yung-Sup  Shin,  American  Bureau  of  Shipping 


Friday,  July  8 

8:00  -  8:30  Coffee  and  doughnuts 

Simulation-Based  Design  Environment 

8:30  -  8:50  Use  of  Reliability  in  Structural  Design 

Mr.  Robert  A.  Sielski,  Marine  Board,  National  Research  Council 

8:50  -  9:25  Virtual  Reality  in  Design  and  Manufacturing 

Prof.  K.-Peter  Beier,  Naval  Architecture  and  Marine  Engineering, 
University  of  Michigan 

9:25  - 10:00  The  Role  of  Simulation  in  Ship  Design:  Some  Cautionary  Examples 

Prof.  Armin  W.  Troesch,  Naval  Architecture  and  Marine  Engineering, 
University  of  Michigan 

10:00-10:15  Break 

10:15  - 10:50  Continued  discussion: 

Dynamic  Loading  Approach  for  Analyzing  the  Ship  Structure 
Dr.  Yung-Sup  Shin,  American  Bureau  of  Shipping 

Workshop  Wrap-up 

1 1 :25  -  12:00  Dr.  Peter  Majumdar,  Office  of  Naval  Research 
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SAFE,  RELIABLE,  AFFORDABLE,  INNOVATIVE  DESIGNS 
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HYDRODYNAMIC  LOADS  TECHNOLOGY 


-  FULL  SCALE  TRIALS 

CPF  (CANADIAN  PATROL  FRIGATE) 

-  MODEL  TESTS 

-  FULL  SCALE  TRIALS  (PLANNED) 
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FRACTION  OF  SHIP  LENGTH 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


0.15  0.2  allen  engleO.25 

IVE  FREQUENCY  (HZ) 

'  r  N&VAI  fiRA  SYffTTUfi  T 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


HYDRODYNAMIC  LOADS  TECHNOLOGY 

PAGE  NUMBER  16 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


2?  CL 
g  CL 

gs 
Q  <  S 

oOO 


12" 


§°5 

Q  IU 

s<s 

P“S 

iu  5  -i 
a2o 

khz 

»-<r 


<5b 

hlZ 
§SW 
z  o  111 
woo 

uj  jy  ° 

nQQ 

oOz 

o°<co 

iu  m  <  s  UJ 
H-  <  CO  X  u. 

O  .  .uiu. 

Z  -r-  oi  I  UJ 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


ALLEN  ENGLE 


HYDRODYNAMIC  LOADS  TECHNOLOGY 


H 

0 

to  5 

w  i 

sc « 

U  o 

c!  Z 

to  ^ 
CD  CD 
to  to 

H 

to  < 

ET  H 
to  CD 

5  D 

25  w 

6  H 

to  HH 

to  z 

O  D 


H 

Z 

W 

S 

CD  ? 

w 

2  > 

2  5 

ri  w 

y  > 

Z  *C 

a  £ 


aj  bj 
w  Q 

a  o 

g  w 

3  S 

O  %  g 

^  §  2 

2wO 

1  g  § 

§  d  S 
s  2  s 

4  P  2 

|gl 

<  g{  CD 

O  H  W 

i-9  CD  Q 
I  I  I 


E 

n  3 

u  « 

p  ^  5 

<  «  3 

3  u  S 

w  53  O 

o  S  o 

Btf  <  gs 

O  Z  to 
£  ><  g 

ffi  Q  < 

.  .  w 

I  I  V 


5 

H 

Z 

I  S 

<  w 
r#s  W  to 

3  >*  to 
5s  «  S 

O  Sz 

H  I  2 

z  g  5 

o  i  2 

g  c  u 
w  z  w 

w  a  h 

Qz  z 

Ph  •“■ 

to.  o  2 

o  a 

hJ  >  53 

tu  w  w 

>  Q  D 


u 

£ 

c 

Q 

ft 

cc; 

O 

£ 

V) 

C3 

£ 

O 

Q 

ft 

cs 

o 

£ 

o 

E 

a 

< 

o 

m 


LARGE-AMPLITUDE  MOTION  AND 
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STRUCTURAL  RESPONSE  PREDICTIONS 
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.  Structural  Design  Approach 


BULK  CARRIERS,  A  CAUSE  FOR  CONCERN' 

Phil  Rynn,  ABS/Houston 


OUTLINE  OF  PRESENTATION 
AND  SPONSORS 


I.  Status  Report  on  IDEAS  Motion  and  Load  System 


A .  System  Development 

1.  ARPA  (1988-90) 

2.  USCG  (1989-94) 
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Figure  2:  Master  Geometry  and  Panel  Distribution  in 
both  Physical  and  Computation  Domains. 
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WAVE  INDUCED  STRUCTURAL  LOADS 
BY  THE  CODE 
SWAN 


by  Paul  D.  Sclavounos 
MIT  Department  of  Ocean  Engineering 


Presentation  to 

ONR  WORKSHOP 

ON  NONLINEAR  SEA  LOADS  AND  SHIP  RESPONSE 
A  BASIS  FOR  SHIP  STRUCRURAL  DESIGN 

College  of  Engineering 
University  of  Michigan,  Ann  Arbor 


July  7&8,  1994 


REVIEW 

OF  SHIP  MOTION 
AND  STRUCTURAL  LOAD 
PREDICTIONS 


BY  THE 

FREQUENCY  DOMAIN  CODE 


SWAN-1 


SWAN2 

transom  hull  mesh 


Figure  10  :  Heave  and  pitch  exciting  forces  and  motion^  of  a  modified  Wigley  model  advancing  at 
Froude  number  F =0.3  through  regular  head  waves. 
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21  TU 


STRUCTURAL  LOADS  EVALUATED  BY 

SWAN 

A  2 


=  F-s)  =  SHEAR  FORCE  1  ABaTT 

5  «  SENDING  MOMENT  frame 


SEAKEEPING  BY  THE 

TIME  DOMAIN  CODE 

SWAN-2 


'5 


distance  travelled=(Ut)/l=0  i 


Figure  5-2:  Transient  wave  elevation  and  body  pressure  contours  for  the  Series 
hull  in  steady  forward  motion  at  T  —  0.2.  started  impulsively  from  rest  at  t  —  0. 


Figure  5-8:  Radiated  wave  pattern  for  the  Series  60  hull  in  forced,  periodic  heave  at 
T  =  0.2  and  encounter  frequency  u >/(g/L)i  =  3.335,  viewed  from  above  and  behind 
the  vessel.  A  snapshot  of  the  steady-state  periodic  wave  pattern  taken  at  the  middle 
of  the  heave  cycle. 


87 


Figure  5-13:  Cross-coupling  added  mass  and  damping  coefficients  for  the  Series  60 
hull  in  heave  and  pitch  at  T  —  0.2. 


Figure  5-14:  Steady  wave  pattern  for  a  transom  hull  at  T  —  0.3,  viewed  obliquely 
from  above  and  behind  the  vessel. 
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Figure  Heave  and  diffraction  wave  patterns  and  linear  body  pressure  patterns 

for  a  transom  hull  at  T  =  0.3  and  encounter  frequency  W(s/lU'=  3.2. 
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7.5 


Figure  6-8:  Magnitude  and  phase  of  the  pitch  response  amplitude  operator  for  the 
Series  60  at  J-  =  0.2. 
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Figure  1:  Description  of  Coordinate  Systems  and  Geometry  of  the  Body 


Figure  2:  Photograph  of  the  Oscillating  Body  and  its  Supporting  Structure 
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Figure  17b;  Comparison  between 
Theory  and  Experiment  ot  the  Verti 
cal  Slam  Coefficient  for  a  Cusped 
body  (FN-3.S371). 
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FULLY  NONLINEAR 
HYDRODYNAMIC  LOADS  USING 
DESINGULARIZED  METHODS 

Robert  F.  Beck 
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Yusong  Cao 
Steve  Scorpio 
Minglun  Wang 


PROBLEM  FORMULATION 


*  Basic  Assumptions: 

1.  Incompressible  and  Invisdd  fluid 

2.  Irrotational  flow 

3.  Surface  tension  neglected 


4  =0  (ISO) 

4-0  (ISO) 


•  • 


FREE  SURFACE 
BOUNDARY  CONDITIONS 


Kinematic  condition: 


(on  F.S.) 

5x 

5t  "Vx 

(on  F.S.) 

8t  y 

(on  F.S.) 

Dynamic  condition: 

f* = -gn  -  \  v*  •  v* + v  •  v*  -  A  -  u0(t)|* 

(on  F.S.) 

where  %- = + v  •  v  is  the  time  derivative 
8t  at 

following 

a  node  moving  with  velocity  v . 


Fixed  Horizontal  Nodes 


and 


||- v<t>.  Vn-U0(t)|n  (orF.S.) 


Material  Nodes  (v = uc(t)i + v<|>) 


and 


^■  =  U0(t)i  +  V$ 


where  xF(t)  =  (xF(t),  yF(t),  zF(t))  is  the  position  vector  of 
a  fluid  particle  on  F.S.  and 


lUA+vov 

Dt  dt 


is  the  material  derivative. 


FULLY  NONLINEAR 
SOLUTION  METHOD 


Time-Stepping  Procedure 


1.  Solve  a  mixed  B VP  at  a  given  instant  of 
time  by  a  desingularized  boundary  integral 
method. 

2.  Integrate  the  nonlinear  free  surface  kinematic 
and  dynamic  conditions  with  respect  to  time. 

3.  For  free  body  problem,  must  compute  g  or 
^  at  present  time  step,  then  integrate  the  body 

at 

equations  of  motion. 


DESINGULARIZED  BOUNDARY 
INTEGRAL  METHOD 


Simple  sources  can  be  used  because  of 
desingularization 

Easy  discretization  (nodes  only,  no  panels) 

Desingularization  distance  related  to  the  local  node 
spacing 

Simple  algorithm,  easy  programming  and  high 
performance  on  supercomputers 

Can  lead  to  0(N)  algorithm 


Figure  1:  Schematic  of  source  and  node  locations 


HYDRODYNAMIC  FORCES 


where 


F>=_JJsPn>ds 

P  3<j>  ,  -  .  .  9<j>  1  rj,  „ , 

r~*  o(t)aTg2~5  *  * 

=-f*~Uotf“8z'iv*'v*+T'v* 

0  =  perturbation  potential  in  moving  coordinate  system. 


APPLICATIONS  OF  THE 
FULLY  NONLINEAR 
DESINGULARIZED  METHOD 


Verification  by  comparison  with  analytic  solutions  for  flows 
generated  by  isolated  singularities  and  bodies  of  simple  geometry 
(Cao,  Schultz,  and  Beck  1991  and  Cao,  Lee,  and  Beck  1992) 

Shallow  water  solitons  due  to  a  moving  disturbance 
(Cao,  Beck,  and  Schultz  1993a) 


Two-dimensional  wave  tank  with  an  adsorbing  beach 
(Beck,  Cao,  and  Lee  1993  and  Cao,  Beck,  and  Schultz  1993b) 


Submerged  spheroid  traveling  at  constant  forward  speed 
(Bertram,  Schultz,  Cao,  and  Beck  1991) 


Two-dimensional  added  mass  and  damping  for  forced  heave,  sway 
and  roll  (  Lee  1992,  and  Beck,  Cao,  and  Lee  1993 ) 

Two-dimensional  heave,  sway  and  roll  motions  of  a  rectangular 
body  due  to  incident  waves  (Cao,  Beck,  and  Schultz  1994) 


Three-dimensional  cylinder  in  forced  heave  (Beck,  Cao,  and  Lee 
1993) 


Exciting  forces  on  a  Tension  Leg  Platform 

Calm  water  resistance  and  added  mass  and  damping  of  a  Wigley 
hull  (Beck,  Cao,  and  Lee  1993,  and  Beck,  Cao,  and  Scorpio  1994) 


Free  Surface  Elevations 
for  a  2-d  tank  with  a  45  degree 
wedge  shaped  body 

with  and  without  spray  damping 
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Freely  Floating  Body  Dynamics 


Hydrodynamic 

problem 


^=f1(xf,4.f,xg,vg) 

^=f2(xf,4.f.xg,vg) 


Euler's  equation 
of  motion 


f“fi-=vG 


dt 


^■=f3[xf,*f,xg,vg,^ 


where  the  state  variables  are  defined  by  the  generalized  vectors: 

XF  =  location  of  free  surface  nodes 
<J>F  =  potential  of  free  surface  nodes 

Xq  =  location  of  vessel  center  of  gravity  in  6-degrees  of  freedom 
Vq  =6  components  of  body  velocity  at  center  of  gravity 


Hydrodynamic  force  at  present  time  step 

Fi=-JJpnids 

s 

ib-¥-u°<-84v^ 

=  -f£"Uo(t)i£-gZ-£V*'V*  +  V'V* 

To  evaluate  ^  or  ~ 

at  5t 

•  Backwards  differencing  for  ^ 

•  Direct  solution  of  ^  BVP 

dt 

v2^J=o  in  fluid 
^=-gn-u0(t)^-iv*  v*-i  onSF 

l(f)  =  ¥(VH'V*)  +  ni|i  0nSH>SB 

Vg)^°  R-»- 


WAVE-INDUCED  FREE  BODY  MOTIONS  IN 
A  TWO-DIMENSIONAL  WAVE  TANK 


Pneumatic  wavemaker  Floating  body  Wave  absorbing  zone 


f.umb') 


f.HmR') 


FJfaH*)  YgfH 


/•,/(«»*)  ,  w» 


M.UmB')  Hr O) 


y(*.*)  = 

where 


Wigley  Hull 


L  =  model  length 

B  =  model  full  beam 

T  =  model  draft 

=  coefficient  for  bow  fullness 
=  0.0,  standard  hull 
=  .2  for  modified  Wigley  hull  m 


For  both  the  standard  hull  and  the  modified  hull  HI,  L/B=10  and  B/T=1.6. 


Desingularization  near  the  leading  edge  of  a  Karman  - 

Trefftz  airfoil 


Desingularized  distance  at  leading  and  trailing  edge  =  0.1  *Ld 
Entrance  half  angle  =  13.5  degrees 


Desingularized  distance  at  leading  and  trailing  edge  =  0.5*Ld 
Entrance  half  angle  =  13.5  degrees 


x 


The  effect  of  desingularized  distance  on  surface  tangential 
velocity  (vt)  for  a  Karman  -  Trefftz  airfoil 


Desingularized  distance  at  the  bow  =>  o.I *Ld 

0.5  *Ld 


Wigley  hull  double  body  solution 
pressure  on  the  forward  half  of  the  body 


Wave  profile  along  the  standard  Wigley  hull 
(Fr  =  0.25,  fixed  sinkage  and  trim) 


wave  amplitude  (m) 


Wave  Elevation  Along  Centerline  and  Wigley  Hull,  Fr«0.25 


x(m) 


0.0005 


•0.0005 

-0.001 


Added  Mass  and  Damping 
in  Heave  and  Pitch 

for 

Modified  Wigley  Hull  III 


F15/(l/2pU02S0) 


F35/(l/2pU02S0) 


Modified  Wigley  hull  HI  -  heave  force  due  to  pitch  excitation 
Fr  =  0.3,  pitch  amplitude  =  1.5\  co/Vg/L  =  2.76642 


FJ5/(l/2pU02S0L) 


Modified  Wigley  hull  HI  -  pitch  force  due  to  pitch  excitation 
Fr  =  0.3,  pitch  amplitude  =  1.5%  co/Vg/L  =  2.76642 


A53/(pVL)  A33/(pV) 
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Comparisons  of  experimental  and  theoretical  added  mass  and  damping 

coefficients  for  forced  heave 
(Wigley  model  ID,  Fr  =  0.3,  z«/L=  0.00833 ) 


A35/(pVL)  A55/(pVL2) 


Comparisons  of  experimental  and  theoretical  added  mass  and  damping 

coefficients  for  forced  pitch 
(Wigley  model  IE,  Fr  =  0.3,  pitch  amplitude  =  1.5°) 


Forced  Oscillation 
of 

Large-Flare  Body 
U0  =  0 


Time  History  of  Force  on  the  Flare  Body  (Numerical) 


(  qO  zd 


t  (sec.) 


Time  History  of  Force  on  the  Flare  Body  (Numerical) 


(•qO  zd 


Time  History  of  Force  on  the  Flare  Body  (Experiment) 


(•qi)  zd 


t  (sec.) 


Time  History  of  Force  on  the  Flare  Body  (Experiment) 


Time  History  of  Force  on  the  Flare  Body  (Numerical) 


Oqi)  m 


Frequency  =  0.597  (Hz) 
Amplitude  =  1.50  (in.) 


Time  History  of  Force  on  the  Flare  Body  (Numerical) 


Frequency  =  0.: 
Amplitude  =  l.t 


LOADS  ASSOCIATED  WITH  THE 
HYDRODYNAMIC  IMPACT  OF 
FLAT  WEDGES  (flat  cylinders) 


William  S.  Vorus 

The  University  of  Michigan 


a)  Zero  Viscosity 

b)  Zero  Compressibility 

c)  Zero  Gravity 

Hydrodynamic  Impact  of: 

1)  Relatively  Flat  Cylinders  (2D) 

2)  of  Otherwise  Arbitrary  (smooth)  Geometry 

3)  Including  Variable  Impact  Velocity 

(with  oscillatory  perturbation) 

4)  Time  Varying  Contour  Shape 

5)  Multiple  Contours 


DEMONSTRATION  IN  TERMS  OF  SELF-SIMILAR  SEMI-INFINITE 
WEDGE  IMPACT 


WEDGE  IMPACT  FLOW 

FIGURE  1 


CYLINDER  PENETRATION  (CW) 

FIGURE  1b 


FREE  CONTOUR  DYNAMIC  BOUNDARY  CONDTHON 


Cp(0  -  0  on  £  £  1 


CAC)  =  l-vfo-v}(0+l‘a 


Lf0 


I 


0  =  l-v/(f)+2za 


b 

I 


ys(c0^0^cvs(o 


This  condition  is  clearly  satisfied  fot  V$(0  =  Vj,  a  constant  in  l£f 
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WEDGE  CONTOUR  KINEMATIC  BOUNDARY  CONDITION 
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Here,  A  =  and  F  is  the  hypergeometric  function  of  one  argument 
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VELOCITY  CONTINUITY 


Remove  singular  terms  as: 
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DISPLACEMENT  CONTINUITY 


SYSTEM  SOLUTION  SUMMARY 
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Continuity  of  Pressure 
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Continuity  of  Velocity 
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Wedge  Vortex  Distribution 
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Wedge  Pressure  Distribution 
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COMPUTATIONS 


Wetting  Factor  ("jet  rise") 
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WEDGE  SOLUTION  CHARACTERISTICS 

TABLE  II 
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JET  THICKNESS  AT  JET-HEAD  TRUNCATION 

FIGURE  8 


Jet  thickness  distribution,  for  Cp  =  0: 


8(s)  =  8. 


f  V.-j  ' 


Ki  zbt 


zbt*s*zit 


For  5(zit)  =  0: 


ADDED  MASS  OF  IMPACTING  WEDGES 
FIGURE  6 
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Wedge  Side-Angle,  p  (deg) 


JET-HEAD  SEPARATION  VELOCITY,  Vs(l.x), 
AND  OFFSET,  bc(x),  versus  Time  ,x 


VS(1,x)/10 
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FUTURE  WORK 


Nonsymmetric  Cylinders 

Linear  Gravity  Flow  Beyond  the  Jet  Head 

Lifting  Surface  (3D)  Corrections 

Extraction  (inverse  impact) 

Rapid  Lateral  Expansion  of  Thin  Cylinders 
(vertical  line-source) 


Gas  Entrapment 


ONR  WORKSHOP 

NONLINEAR  SEA  LOADS  AND  SHIP  RESPONSE:  A  BASIS 

FOR  SHIP  STRUCTURAL  DESIGN 
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The  University  of  Michigan 


Issues  Related  to  the 

Hydrodynamics  of  Planing  Hull  Seakeeping 


The  University  of  Michigan 
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The  University  of  Michigan 


The  University  of  Michigan 


The  University  of  Michigan 


Validation  of  numerical  model 


The  University  of  Michigan 


Deadrise:  20  degrees;  Trim:  4  degrees; 
Mean  wetted  length/beam  ratio:  2.56. 


Validation  of  numerical  model 


The  University  of  Michigan 


MOTION  DEPENDENT  WETTED  LENGTH 


The  University  of  Michigan 


Evaluation  hydrodynamic  loads 


The  University  of  Michigan 
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Moment  stiffness  due  to  pitch 


The  University  of  Michigan 


Hydrodynamic  inertial  force  due  to  heave 

(First  order) 


The  University  of  Michigan 


Hydrodynamic  damping  force  due  to  heave 


The  University  of  Michigan 


namic  inertial  moment  due  to  heave 


The  University  of  Michigan 


Hydrodynamic  damping  moment  due  to  heave 

(First  order) 


The  University  of  Michigan 


Hydrodynamic  inertial  force  due  to  pitch 

(First  order) 


The  University  of  Michigan 


SHIP  STRUCTURES 
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FIGURE  70  LIFETIME  LOADS  ANALYSIS  METHODOLOGY 


OPERATIONAL  PROFILES 


vo  m  m  ^  © 
^  A  ^ 

vo 


FIGURE  28  TYPICAL  FILTERED  TIME  HISTORIES  OF  WAVE  INDUCED  BENDING, 
SHOWING  ORDINARY  WAVE  AND  WHIPPING 
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FIGURE  38  TYPICAL  WEI  BULL  PLOTS  OF  CG-61  FULL  SCALE  MIDSHIP  WHIPPING  MOMENTS 
C)  VERTICAL  BENDING.  8EA  8TATE  8. 20  KNOTS.  HEAD  SEAS 


LN(WHIP  MOMENT/CHARACTERISTIC  VALUE) 


FIGURE  88  TYPICAL  WEIBULL  PLOTS  OF  CO-61  FULL  SCALE  MIDSHIP  WHIPPING  MOMENTS 
d)  VERTICAL  BENDING.  SEA  STATE  6, 25  KNOTS.  HEAD  SEAS 


BASED  ON  HEAD  SEAS 
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FIGURE  57  CG-47  CLASS  MODEL  BASED  VERTICAL  BENDING  MOMENT  ROOT  RAO'S  AS  A 

FUNCTION  OF  SEA  STATE 


b)  CG-47  CLASS  PHASE  ANGLE  -  BON  FLAKE  IMPACTS 

FIGURE  27  TYPICAL  BENDING  MOMENT  TIME* HISTORIES  SHOWING  PHASE 
ANGLES  BETWEEN  ORDINARY  WAVE  AND  WHIPPING 
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INTEGRATED  SHIP  STRUCTURAL  DESIGN 

METHODOLOGY 


Prof.  Owen  Hughes 

Aerospace  &  Ocean  Engineering  Department 
Virginia  Polytechnic  Institute  and  State  University 

Tobin  R.  McNatt 
Proteus  Engineering 
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Program  Objectives 

•  Leveraging  of  past  work  in  ship  structural  design 
conducted  in  the  US  into  a  practical,  computer-based 
system. 

.  Resolution  of  important  technical  requirements  related  to 
fatigue,  ultimate  strength  and  reliability  which  will 
facilitate  the  implementation  of  a  reliability-based  design 
system. 

•  Utilization  of  ongoing  efforts  in  predicting  ship  motions 
and  loads. 

•  Preserving  the  technical  leadership  of  the  US  in  ship 
structural  design  capabilities. 

.  Providing  advances  in  the  safety  and  cost  effectiveness  of 
military  and  commercial  ship  structures,  offering  dual-use 
of  the  integrated,  computer-based  structural  design  tool 
which  will  be  produced  by  the  research  efforts. 

.  Enabling  the  rational,  reliability-based  structural  design  of 
ships  from  first  principles,  so  that  new  or  advanced 
geometries  (e.g.  double  hulls,  SWATHs)  can  be  practically 
designed  by  non  research  groups  such  as  shipyards  and 
naval  architects. 
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Overview  of  Research  and  Development  Program 
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Q  Ship  Structural  Design  Technology  Development  ^ 
Structural!  Disciplines  Focus  Group 
(Task  1.  Fatigue  Design 
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Task  2.  Ultimate  Strength 


|(Task3.  Damage  Tolerance 
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(Task  4.  Slam-Induced  Structural  Response 
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Structural  Design  System  Focus  Group 


(Task  5  Structural  Reliability  Technology 


(  Tatfc  ».  Interfacing  Hydro  -  Structures  Codes 


(^  Tit*  I.  Multi-Disciplinary  Design  Optimization 


(  Task  7.  Integrating  Structural  Design  System 
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Year  4 


Year  5 


Virginia 

|J]Tfech 

Years  One/Two  Tasking 

Task  1.  Develop  a  Fatigue  Design  Procedure 
Task  2.  Improve  Ultimate  Strength  Analysis  Methods 
Task  2.  Develop  Damage  Tolerance  Design  Technology 
Task  4.  Computation  of  Slamming-lnduced  Structural 
Task  5.  Structural  Reliability  Technology 


Virginia  , 

mifech 

Task  1.  Develop  a  Fatigue  Design  Procedure 

•  Reformulate/extend  the  analysis  method  and  the  design 
method  to  include  fatigue 

•  Define  fatigue-related  load  data  as  needed  from  a  Hydro 
program  (e.g.  the  influence  of  whipping) 

•  Define  a  modular  interface  with  such  a  program 

.  Extend  the  stress  analysis  to  get  cyclic  stress  transfer 
functions;  initially  for  stiffeners  and  frames 


Virginia 

|||Tech 

Task  2.  Improve  Ultimate  Strength  Analysis  Methods 

.  Member  Level  -  top  priority  need:  a  better  model  for 
flexural-torsional  buckling  of  beams 

•  Overall  Level  -  further  improve  the  current  analysis  model 
to  better  account  for  post-buckling  stiffness 


Virginia 
Ijg!  Tech 


Task  3  Develop  Damaae  Tolerance  Desian  Technolo 


•  Define  relevant  load  conditions  for  damage  tolerant  design 
(e.g.  grounded  and/or  flooded  conditions) 

•  Develop  detailed  plan  for  incorporation  of  Damage 
Tolerance  Approach  into  design  practice 


Virginia  , 
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Task  4  Computation  of  Slamminq-lnduced  Structural 
Responses 

•  Define  local  and  global  slamming-related  load  data  needed 
from  Hydro  program(s) 

•  Later  phases  will  address  the  computation  of  structural 
response  to  slamming  and  the  interfacing  between  slam 
load  prediction  codes  and  structural  response  codes 


Virginia 
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Task  5  Structural  Reliability  Technology 

•  Define  the  basic  technology  requirements  for  a  reliability- 
based  format  for  ship  structures 

•  Assess  the  state-of-art  for  the  technology  requirements 
and  determine/prioritize  the  development  needs 

•  Develop  a  plan  for  ensuing  years'  efforts 

•  Identify  and  establish  liaison  with  current  relevant 
research 
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Implementation  of  Results 


•  Overall  Objective:  To  incorporate  these  developments, 
together  with  the  other  components  of  design,  into  an 
integrated  structural  design  system  which  can  be  used  by 
the  ship  design  community 

•  ONR  has  decided  that  MAESTRO  will  be  the  vehicle  for 
implementing  the  results  into  an  existing  structural  design 
tool 

•  MAESTRO  integration  with  other  design  tools  has  been 
selected  under  the  Maritech  Modular  Tanker  Consortium 
project 


Our  team  will  also  be  participating  in  the  US-Norway 
research  project  Dynamic  Analysis  of  Surface  Ships. 
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Structural  Load  Cases 


Dynamic  Loading  Approach 


Fatigue 

Assessment  n — "t - >  Structure 


Introduction  of  DLA 

For  Added  Margin  of  Safety 
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Dynamic  Loading  Approach  (DLA) 


Equivalent  Wave  in  Partial  Loading  (67%  LOAD(C)) 
DLP  :  Maximum  Hogging  Condition 


THE  DYNAMIC  LOADING  APPROACH  (DLA) 

FOR 

ANALYZING  SHIP  STRUCTURE 


Presented  at  ONR  Workshop  on 
Nonlinear  Sea  Loads  and  Ship  Response 
Ann  Arbor,  Michigan 

7-8  July  1994 


Yung  S.  Shin 


AMERICAN  BUREAU  OF  SHIPPING 
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A  Typical  3-D  FEM  Model  of  a  Tanker 
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3-D  Global  FEM  Analysis 
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Loading  for  FEM  Model 


Internal  Tank  Pressure  Distribution 


Pressure  Components  due  to 

•Vapor  Pressure 

•Roll  and  Pitch  Inclination 

•Accelerations 


External  Pressure  and  Distribution 


Pressure  Components  due  to  | 

•  Wave  I 

•Vertical  Motion  I 

•Lateral  Motion 

•Roll  I 


ROLL  ANGLEtDEG. 


N.  BALLAST 
33%  PARTIAL 


ROLL  ANGLE  V. S.  SHIP’S  DISPLACEMENT 
(AT  MAX.  ROLL  CONDITION) 


Long  Term  Response  and 
Equivalent  Wave  System 


Short  Term  Response 


WAVE  SPECTRUM 
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FRF  OF  DESIGN  LOAD  PARAMETER 
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RESPONSE  SPECTRUM  OF  DIP 


Representative  Frequency  Response  Function 

_  Vertical  Bending  Moment 


Representative  Frequency  Response  Function 

Roll  Motion 


Equations  of  Motion 

6-Degrees  of  Freedom  Motion 


•  Surge 

•  Sway 

•  Heave 

•  Pitch 
•Roll 
•Yaw 


[{Mfit  +  Ajk)i}k  +  Bjkik  +  CjkVk] 

=  FJ«;  j  =  1...6 


fji  =  surge 
i)2  =  sway 


773  =  heave 
ru  =  roll 


iji  =  pitch 
v  6  =  yaw 


Critical  Load  Cases 


Load  Components 


Representative  Locations  of  Acceleration 


APPLICATIONS  AND  BENEFITS 


VIRTUAL  PROTOTYPING : 

-  reduce  or  eliminate  the  need  for  costly  and  time 
consuming  physical  prototypes 

-  rapid  virtual  prototyping  allows  for  fast  creation  and 
modification  of  prototype 

ENGINEERING  ANALYSIS : 

Integration  of  analysis  results  (e.g.,  CFD,  FEM)  with 

virtual  prototypes 

OPERATIONAL  SIMULATIONS : 

-  direct  involvement  of  humans  for  ergonomic,  human 
factors,  and  performance  studies 

-  simulation  of  assembly,  production,  and  maintenance 
tasks  reveal  problems  at  an  early  stage  of  the 
design  process 

-  training  for  operation,  maintenance,  safety 

CONCURRENT  ENGINEERING  : 

-  explore  all  aspects  of  a  design  or  a  process  by  an 
Interdisciplinary  engineering  team 

-  shared  virtual  environments  allow  for  participation 
from  remote  locations 

-  VR  as  an  Integrating  tool  for :  design  -  engineering  - 
analysis  -  production  planning  -  manufacturing  - 
marketing  &  sales  -  maintenance  •  training 

OVERALL  BENEFITS: 

savings  In  cost  -  savings  In  time  -  reduced  design 

cycle  -  Improved  market  response  -  better  product 


THE  IMMERSIVE  EXPERIENCE  OF  VR 


•  Convincing  Illusion  of  Being  Fully  Immersed  in  an 
Artificial  Three-Dimensional  World 

•  Depth  Perception  through  Stereo  Viewing 

•  Full  Look-Around  &  Walk- Around  Capability 

•  Full  Scale  Representation  of  Virtual  World 

•  Realistic  Interactions  with  Virtual  Objects 

•  Strong  Sense  of  Realism  and  Spatial  Perception 


ASSETS  FOR  DESIGN  AND  MANUFACTURING 


•  Optimal  Analysis  Tool  for  Spatial  Problems  Involving 
Complex  Three-Dimensional  Geometry 
(arrangements,  mechanical  systems,  abstract  systems) 

•  Realistic  Integration  of  Humans  with  Virtual  World 
(especially  effective  if  humans  are  part  of  a  system) 

•  Optimal  Communication  &  Demonstration  Tool 


Figure  1 .  A  Head-mounted  display  (HMD)  provides  fully  immersive  virtual  reality. 


Figure  2.  Data  glove  with  flexion  sensors  and  tracking  receiver. 


THE  ENABLING  TECHNOLOGIES  OF  VR 


•  Head-Mounted  Display  (HMD) 

•  Motion  Tracking  System 

•  Image  Generation  System 

•  Interactive  Input  Devices  (Gloves,  Suits) 

•  Tactile  and  Forced  Feedback 

•  Additional  Component  Technologies 

Eye  Tracking  Systems 
Telepresence  Technologies 
Directional  3D  Sound 
Voice  Recognition  /  Speech  Synthesis 

•  Alternative  Display  Technologies 

Head-Coupled  Display  (HCD) 

See-Through  HMD/HCD 
Shutter  Glasses 
Large  Screen  Projection 
Retinal  Laser  Display 


Klaus-Peter  Beler 


University  of  Michigan 


VIRTUAL  REALITY  (VR) 
VIRTUAL  ENVIRONMENTS  (VE) 

IN  DESIGN  AND  MANUFACTURING 


•  The  Enabling  Technologies  of  Immersive  VR 

•  Applications  In  Design  and  Manufacturing 

•  Virtual  Prototyping  (case  study:  automotive  Interiors) 

•  Augmented  Reality  In  Assembly  and  Maintenance 

•  The  University  of  Michigan  Virtual  Reality  Laboratory 


CMS-SSC  STRUCTURAL 
RELIABILITY  THRUST 

OTHER  RELIABILITY  PROJECTS 
SSC  363  -  1992 

Uncertainties  in  Estimating  Loads  and  Load  Effects  on 

Marine  Structures 
Estratos  Nikoladis 

Developed  estimates  of  bias  and  uncertainty  in  loads  and 

load  effects 

SSC  371  -  1993 

Establishment  of  a  Uniform  Format  for  Data  Reporting  of 
Structural  Material  Properties  for  Reliability  Analysis 

Fleet  Technology  Limited 

Developed  a  standard  format  so  that  the  properties  of 
materials  will  be  available  in  a  probabilistic  format 

SR  -  1338 

Uncertainty  in  Strength  Models  for  Marine  Structures 

Owen  Hughes 

Objective  -  Quantify  bias  and  uncertainty  in  structural 
strength  formulations  in  order  to  evaluate  safety  margins 

and  derive  design  criteria. 

SR-1344 

Assessment  of  Reliability  of  Existing  Ship  Structures 

Alaa  Mansour 

Will  estimate  the  reliability  levels  associated  with 
important  failure  modes  for  several  existing  ships 


CMS-SSC  STRUCTURAL 
RELIABILITY  THRUST 


PHASE  m  -  SR-1345 

Probability-based  Ship  Design:  Implementation  of  Design 

Guidelines  for  Ships 
Alaa  Mansour 

Project  began  in  May  1994  with  the  objective  of  developing 
ship  structural  design  procedures  that  are  reliability-based. 

PHASE  IV  -  SR-1362  (New  Phase) 
Probability-based  Design,  Synthesis  of  the  Reliability 

Thrust  Area 

Project  to  begin  in  FY  1994 

Have  a  group  of  experts  produce  a  summary  of  the  state  of 

the  art 

PHASE  V  (Old  Phase  IV) 
Probability-based  Design:  Novel  Hull  Forms  and 

Environments 

Will  extend  previous  work  to  unconventional  hull  forms 
and  to  unusual  load  situations 


CMS-SSC  STRUCTURAL 
RELIABILITY  THRUST 

1983  DESIGN,  INSPECTION  AND  REDUNDANCY 

SYMPOSIUM 

Recommended  a  program  of  several  projects  to  determine 
and  unify  reliability  of  marine  structural  systems 

June  1987  ad  hoc  Reliability  Committee 
Developed  Four-Phase  Reliability  Thrust 

SSC-351  -  1990 

Alaa  E.  Mansour  and  Paul  F.  Wirsching 
A  tutorial  on  structural  reliability  theory 
The  basis  for  SSC  reliability  work 

PHASE  I  -  SSC  368  -  1993 
Probability  Based  Ship  Design  Procedures  - 
A  Demonstration 
Alaa  Mansour 

Demonstrated  the  use  of  probability-based  design  in  the 
analysis  of  a  tanker  and  a  combatant  ship 

PHASE  II  -  SSC  373  -  1994 
Probability-based  Ship  Design:  Loads  and  Load 

Combinations 
Alaa  Mansour 

Developed  standard  loads  necessary  for  a  probability-based 

design 


COMMITTEE  ON  MARINE 
STRUCTURES 
AND 

SHIP  STRUCTURE  COMMITTEE 
STRUCTURAL  RELIABILITY  THRUST 


Robert  A.  Sielski 
Marine  Board 
National  Research  Council 

July  7,  1994 
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Stronger,  More  Robust,  and  Longer  Life  Vessels 


Acceptance  Criteria 


►  20-year  Life 


Design  Procedure  Using  Dynamic  Load  Approach 


KAXIMIN  IXt,  StBO,  I-M  STRESS  (IG/CX**2> 


A  Typical  Fine  Mesh  FEM  Model  of  a  Tanker 
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Fine  Mesh  FEM  Analysis  of  Local  Structure 


VIRTUAL  PROTOTYPING 


Replace  costly  and  time-consuming  physical  mockups 

Rapid  Virtual  Prototyping:  fast  creation  and  modification 

Create  virtual  prototype  from  existing  CAD/CAM  data 

Apply  rendering  algorithms,  lighting  models,  texture 

mapping,  and  other  techniques  for  realistic  appearance 

Realistic  Interactions  with  prototype  via  data  glove,  etc. 

Combine  virtual  display  with  physical  elements  if  correct 

forced  feedback  Is  needed  (e.g.,  simplified  seating  buck) 

Examples  of  Extended  Functionality  : 

-  employ  transparent  display  techniques  for  Inspection 
of  hidden  components 

-  use  prototype  already  for  the  analysis  of  Incomplete 
designs  (e.g.,  with  parts  floating  In  space) 

-  superimpose  design  alternatives  for  comparison 

-  allow  for  Interactive  design  modifications  with 
immediate  feedback 

Usage  of  Virtual  Prototypes  : 

-  Design  analysis  (clearances,  packaging  efficiency, 
connectivity,  motion  characteristics,  collision, ...) 

-  Human  factors  studies  (visibility,  reachability, 
accessibility,  comfort,  human  performance,  appeal, ...) 

-  Base  for  other  VR  applications  (engineering  analysis, 
operational  simulations,  concurrent  engineering, 
shared  virtual  environments,  training,  marketing, ..) 


CREATION  OF  VIRTUAL  PROTOTYPES 


GEOMETRY  DEFINITION  : 

CAD/CAM  Model :  mathematical  description  through 
surface  modeling,  solid  modeling,  and  related  methods 

Virtual  Model :  approximation  by  computer  graphics 
primitives  (mainly  polygons  and  polygon-meshes) 


PROCESS : 

Access  CAD/CAM  database  and  extract  geometry 

Approximate  boundaries  by  polygons  (tessellation) 

Simplify  geometry  :  reduce  number  of  polygons  for 
graphics  performance,  decide  on  level  of  detail 

Edit  geometry  :  Identify  and  remove  unnecessary 
geometry,  correct  faulty  geometry,  Incorporate 
textures 

Define  display  characteristics  :  color,  reflection 
characteristics  (material  properties),  transparency, 
lighting  configuration,  rendering  method, .... 

Model  the  behavior  of  prototype  (define  Interactive 
feedback  mechanisms  and  trigger  events,  dynamic  and 
other  responses,  motion  restrictions,  etc.) 

Calibrate  virtual  environment  with  user,  data  glove  and 
physical  elements 


73  VW  Super  Beetle-H 


10364  Vertices 
10514  Polygons 


Polygonal  representations 
with  different  levels  of  accuracy 


3 


Figure  10.  Polygonal  approximation  of  a  panel  derived  from  a  CAD/CAM  model  (top) 

and  simplified  VR  geometry  (bottom). 


File  Edit  Selection  Editors  kUnips  Uqhts  Options 


The  process  of  generating 
a  uirtual  prototype  from  a  CRD/CRM  model 


VIRTUAL  REALITY  LABORATORY 


THE  UNIVERSITY  OF  MICHIGAN  -  COLLEGE  OF  ENGINEERING 
Department  of  Naval  Architecture  and  Marine  Engineering 


Created:  April  1993  through  an  initial  grant  from  Chrysler 
with  cost  sharing  by  The  University  of  Michigan 

Director:  Klaus-Peter  Beier 


Focus:  VR  as  an  integrating  concurrent  engineering  tool  for : 

design  -  engineering  -  analysis  -  production  planning 
manufacturing  -  marketing  &  sales  -  maintenance 


Goals: 


Advance  VR  technology  (cross  the  threshold  of  usability) 
Develop  methods,  algorithms,  and  software  concepts 
Prove  usefulness  through  demonstration  projects 
Assist  industry  with  the  introduction  of  VR 


Laboratory  Equipment : 

•  SGI  Onyx  with  2  processors  and  VTX  graphics 

(currently  upgrading  to  4  processors  and  Reality  Engine  graphics) 

•  Boom2C  from  Fakespace  (currently  upgrading  to  Color  Boom3C) 

-  DataGlove  from  VPL,  Isotrack  from  Polhemus 

•  IBM-RS/6000  with  CAD/CAM  system  CATIA 

•  SGI  lndigo2/Extreme,  SGI  Iris,  SGi  Indy,  HP  and  Sun  workstations, 
Macintoshes,  Scanner,  Video  production  equipment,  others, .... 


VIRTUAL  REALITY  LABORATORY 


ONGOING  PROJECTS : 

«  Virtual  Prototyping  of  Automotive  Interiors 

-  for  design  analysis  and  human  factors  studies 

-  sponsored  by  Chrysler  Corporation 

•  Virtual  Environments  as  an  Analysis  Tool  In  Computational 
Fluid  Dynamics  and  Crash  Simulations 

-  focus  on  high  performance  computing  applications 

-  sponsored  by  DoE  (CRADA) 

-  In  cooperation  with  five  National  Laboratories  and  with 
Chrysler,  Ford,  and  General  Motors 

•  Virtual  Interior  Arrangements  for  Sailing  and  Motor  Yachts 

Department  of  Naval  Architecture  and  Marine  Engineering 

•  Virtual  Model:  Integrated  Technology  Instruction  Center 

future  home  of  the  Virtual  Reality  Laboratory,  under  construction 


IN  PREPARATION  : 

•  Augmented  Reality  in  Simulation-Based  Design 

-  applications  In  design,  assembly,  and  maintenance 
•  sponsored  by  ARPA 

-  to  be  Integrated  with  Chrysler  project  on  virtual  interiors 

•  Shared  Virtual  Environments 

-  simultaneous  immersion  of  users  on  both  sides  of  the  Atlantic 

-  in  cooperation  with  the  Computer  Graphics  Center  (ZGDV) 
Darmstadt,  Germany,  Ford/US,  and  Ford/Europe 

•  Virtual  Diving 

-  virtual  models  of  underwater  terrain,  structures,  ship  wrecks 

-  training  of  divers  and  ROV  operators  (M-ROVER) 

-  planning  of  underwater  operations 


The  University  of  Michigan 
Ann  Arbor,  Michigan 


THE  ROLE  OF  SIMULATION  IN  DESIGN 


RATIONALE: 

Identification  of  design  parameters 
that  are  critical  to  performance. 


The  University  of  Michigan 


dtAUdif  jOifjjfcf 


(a)  Single  branch  (b)  Multiple  branches 


The  University  of  Michigan 


The  University  of  Michigan 


The  University  of  Michigan 


SIMULATED  FORCED  MOTIONS 


Background  on  Damaged 
Stability 


The  University  of  Michigan 


A  Single  Degree  of  Freedom  Model  for 

Vessel  Capsizing 
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The  University  of  Michigan 


Integrity  Plots  for  Short  Time  Exposure  With 
Increasing  Wave  Height 


Distribution  of  State  Variables  in  Phase  Space 
Leading  to  Failure  Based  upon  Simulation 


The  University  of  Michigan 


Prediction  of  Failure  Based  Upon 
Phase  Flux  Concepts _ 
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The  University  of  Michigan 


Summary  and  Conclusions 


incorporate  multiple  degrees  of  freedom, 
large  bias  and  random  excitation. 


